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1 Introduction 



The rise of quantum estimation theory dates back to the late 1960s - 1970s 
(see [8], [9]). It was initially developed as an adequate mathematical frame- 
work for design of optimal receiver in quantum communication channels, and 
later turned out to be relevant also for clarification of some foundational is- 
sues of quantum measurement. New interest to quantum estimation theory 
was brought by the development of high precision experiments, in which re- 
searchers operate with elementary quantum systems. In such experiments 
quite important is the issue of extracting the maximum possible information 
from the state of a given quantum system. For example, in currently dis- 
cussed proposals for quantum computations, the information is written into 
states of elementary quantum cells - qubits, and is read off via quantum mea- 
surements. From a statistical viewpoint, measurement gives an estimate for 
the quantum state, either as a whole, or for some of its components (param- 
eters). The main concrete models of interest in quantum estimation theory 
considered up to now fall within one of the following classes: 

1) Parametric models with a group of symmetries [9], [12]. In particular, 
the models with the shift or rotation parameter are strictly relevant to the 
issue of canonical conjugacy and nonstandard uncertainty relations, such as 
time-energy, phase-number of quanta, etc. It is these models which will be 
our main concern in this paper. 
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2) The full model, in which the multidimensional parameter is the quan- 
tum state itself, i. e. we are interested in estimation of completely unknown 
quantum state. Although in finite dimensions it is a parametric model with a 
specific group of symmetries, it deserves to be singled out both because of its 
importance for physics and of its mathematical features. Especially interest- 
ing and mostly studied is the case of the qubit state, with the 3-dimensional 
parameter varying inside the Bloch sphere. Asymptotic estimation theory 
for the full model in the pure state case was developed in [5], [7], and for 
mixed states in [14]. 

On the other hand, the full model, especially in infinite dimensions, be- 
longs rather to nonparametric quantum mathematical statistics, which is at 
present also in a stage of development. In this connection we would like 
to mention the method of homodyne tomography of a density operator in 
quantum optics [3], [6]. 

3) Estimation of the mean value of quantum Gaussian states. This is 
a quantum analog of the classical linear "signal+noise" problems, however 
with the noise having quantum-mechanical origin. This model was treated 
in detail e.g. in [9]. 

An important distinctive feature of quantum estimation appears in con- 
sideration of series of independent identical trials of a quantum system. In 
the paper [10], devoted to asymptotics of estimation of a shift parameter in a 
quantum state, it was observed that entangled covariant estimates in models 
with independent multiple observations can be more efficient than unentan- 
gled ones. For arbitrary locally unbiased observables this was demonstrated 
on the full model in [15], [5], [7], [4]. This property is a statistical counterpart 
of the strict superadditivity of Shannon information for quantum memoryless 
channels (see. §5.1 in [12]). 

In the present paper we develop further the asymptotic theory of esti- 
mation of a shift parameter in a quantum state to demonstrate the relation 
between entangled and unentangled covariant estimates in the analytically 
most transparent way. After recollecting basics of estimation of shift param- 
eter in Sec. 2, we study the structure of the optimal covariant estimate in 
Sec. 3, showing how entanglement comes into play for several independent 
trials. In Sees. 4,5 we give the asymptotics of the performance of the optimal 
covariant estimate comparing it with the "semiclassical" unentangled covari- 
ant estimation in the regular case of finite variance of the generator of the 
shift group. Sec. 6 is devoted to estimation in the case where the regularity 
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assumption is violated. 

In this paper we deal with the case of pure states. It is in fact yet another 
distinctive feature of quantum estimation that the complexity of the problem 
increases sharply with transition from pure to mixed states estimation. In 
fact, estimation theory for mixed states is an important field to a great extent 
still open for investigation. Another simplifying factor in our model is that 
the parameter is one-dimensional. It is well known that in the quantum case 
estimation of multidimensional parameter involves additional problems due 
to the non-commutativity of the algebra of quantum observables, see e. g. 
[9]. 

2 Covariant Estimates 

Let Ti be a Hilbert space of observed quantum system. The states of the 
system are described by density operators S in H, and observables with 
values in a measurable space X - by resolutions of the identity or probability 
operator valued measures M on X (see, e. g. [9], [12]). Let G be either the 
real line R (the case of a displacement parameter) or the unit circle T (the 
case of a rotational parameter), and let x — > V x — e~ zxA ; x G G, be a unitary 
representation of the one dimensional Abelian group G in H. The spectrum 
A of the operator A is contained in the dual group G, which is identified 
with K. in the case G = R and with the group of integers Z when G — T . 
We consider the problem of estimation of the shift parameter G G in the 
family of states 

S g = e- WA Se t6A , 

where S is the inital state, assumed to be known. By estimate we call arbi- 
trary observable with values in G. 
The estimate M is covariant if 

V*M(B)V X = M(B -x) for B C G, x G G, (1) 

where B — x := {y; y + x G B}, and in the case G = T addition is modulo 2n. 

A necessary and sufficient condition for existence of a covariant observable 
is that the spectrum of A is absolutely continuous with respect to the Haar 
measure in the dual group G [11], which we assume from now on. In the case 
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G — R this means that A has Lebesgue spectrum, while in the case G = T 
(with G — Z) this poses no restrictions. We introduce the operators 

U y := [ e iyx M(dx); yea (2) 
Jg 

Then (1) reduces to the Weyl relation 

U y V x = e ixy V x U y ; xeG,yeG, (3) 

in which, however, the operator U y , is in general non-unitary. In this sense 
the observable M is canonically conjugate to the observable A. 
Introducing the characteristic function 

ipf(X) = E S U X = [ e iXx ^(dx), XeG, 
Jg 

one has the following uncertainty relation for the generalized canonical pair 
(A : M) [10] 

Af(X) ■ D S (A) > \ (4) 

where Ag f (A) := \~ 2 (\)\~ 2 — 1), A ^ 0, is a functional measure of 
uncertainty of the covariant observable M in the state S. If G — R and M 
has finite variance D S (M), then lim A ^ Af f (A) = D S (M), so that from (4) 
follows the generalization of the Heisenberg uncertainty relation 

D S (M) ■ D S (A) > l -. (5) 

Of the main interest are covariant estimates having the minimal uncer- 
tainty. The following result describes them in the case of a pure initial state 
S = where ip is a unit state vector in H. Then 

S fl = e-*»Me* A . 

Since the spectrum of A is absolutely continuous, we have the direct integral 
spectral decomposition 

H= I ®H{\)d\ (6) 

J A 
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diagonalizing the unitary group { e t9A } . This means that for any vector 



\<f) 



®\tp{\))d\, \<p(\)) eH(\), 



one has 



J8A 



M = L 



e idX \<p(\))d\, 6 G G. 



We agree to extend y?(A) to the whole of G by letting it zero outside of A. 
Theorem 1 [10]. For arbitrary covariant estimate M 



<¥*(A) := / + ^ II d\', 



G 



Af(X) > A* (A) := A- 2 (v^(A)- 2 -l). 



In the case G 



D S (M) >D^: 



(7) 
(8) 

(9) 



provided the right hand side is defined and finite. 

The equalities are attained on the optimal covariant observable M* which 
is given by the following kernel in the direct integral decomposition (6) 



M*(dx) = 
Notice that the operator 

P* := / M*(dx) 

JG 



(r-A) \m) ma')i 



ll^(A)||||V(A')|| 



dx 
2?r' 



□ 



(10) 



5(\' - A) 



l^(A))(^(AQ| 

\mx)\\um\ 



(ii) 



is a projection onto the invariant subspace of the group {e* 6 ^} generated by 
the vectors e i6A \ip ),0 e G. Indeed, P*e ieA \ip ) =e ieA \ip) and if the vector \<p ) 
is orthogonal to all of these vectors, then f A e tXe {<p(X)\ip(X) )d\ = 0, 6 G G, 
and ((p(X)\ip(X)) = for A G A, hence P*\<p) = 0. Thus M* is in general 
subnormalized, and one has to extend it to the orthogonal complement of 
= P*TL to obtain an observable. Independently of the extension, it has the 
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following probability density in the state S$, given by the Fourier transform 
of the characteristic function (7 ), 



p* e (x) 



1 

2^ 



U{\)\\d\ 



In [11] it is shown also that M* minimizes the average deviation 1Z(M) = 
fW(x — 9)p^f(x)dx where W is an arbitrary continuous conditionally neg- 
ative definite function. 

Similarly to (11), we obtain for future use 



U, 



*y 



M*(dx) 



G 



S(y + X'-X) 



U(X)\\UM\\ 



(12) 



Example [8], [9]. Consider the case H =L 2 (R), where A acts as multipli- 
cation by the independent variable A (momentum representation). Then 6 is 
the position displacement parameter. In this case T~t(X) — R, ip{X)/ ||^(A)|| = 
e ia(A) an( j ^q-j j g an or thogonal resolution of the identity in TC . This is most 
easily seen by verifying that the operators (12) form a unitary group. This 
resolution of the identity is the spectral measure of the selfadjoint operator 



/l d 
XM*(dx) = - — \y =0 U*y. 



(13) 



According to (12) the action of this operator on the vector (p is given by 



where Q 



dX 



dX 

is the position operator in the momentum representation. 
Thus the optimal covariant observable is position observable up to the gauge 
transformation compensating the phase of the state vector \ip) in the mo- 
mentum representation. In case the argument a (A) is absolutely continuous 
one has further 

Q*y?(A) = [Q + a'(X)]<p(X) 
on an appropriate domain of functions <f(X). 
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3 Multiple Observations 

Now we consider the problem of estimation of the shift parameter 9 G K in 
the family of states 

Sf n = S e © . . . © S e (14) 

in the Hilbert space 'H® n = 7i <S> ■ ■ ■ <S> Ti (n-fold tensor product which 
corresponds to n independent observations). Here Sg = e~ ldA Se z9A , where 
S = \ip)(ip\ is a pure state. 

The family (14) is covariant with respect to the unitary representation of 
the group of shifts of R 

9 -> exp{-i9A (n) ), A {n) = A©... ©/ + ... + /©... ©A 

in 1-t® n . The corresponding direct integral decomposition reads 

where 

H (n) (A) = J © [W(Ai) © ... © H{K-i) © W(A - Ai - . . . - A n _0] dAi . . . dA n _L 

For symmetry of notations we shall also denote this as 

/ © [W(Ai) © . . . © H(\ n )} tf-V, 
Jr(A) 

where 

r(A) = {(A 1 ,...,A n ):A 1 + ... + A n = A}. 
Then |^® n ) = / ©|V> (n) (A))<U, where 

|^ (n) (A)) = / © [|V(Ai)> © ... © |V(An)>] tf-V, 
Jr(A) 

so that 

||^»)(A)|| 2 = / © [|l^(Ai)|| 2 © ... © ||^(A„)|| 2 ] cT-V. (15) 
Jr(A) 
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With these relations in mind, the optimal covariant observable is given by 
the formula (10) 



Mi n \dx) 



= e ix(X'-X) 



|^W(A)) (^(AQ | 1 dx 
||^)(A)||||^)(A0||J 2^ 



(16) 



However in this case the projection 




Mi n \dx)= (J(A'-A) 



\^ n \\)){^(\')\ ' 
||^)(A)||||^)(A0IL ' 



n > 1, 



(17) 



cannot be equal to the identity operator: in any case it projects onto a sub- 
space lying in the subspace of vector functions symmetrically depending 
on Ai, . . . , A n . 

Let us denote p(\) = ^(A)!) 2 and p( n \\) = ^^(A)!! the probabil- 
ity densities of the observables A and in the states S and S® n re- 
spectively. Then (15) implies that p^ n \X) is the n— th convolution of p(A), 
which we denote as p^ n \\) = p(X)* n . We assume that p(X) and hence p^ n \\) 
are differentiable functions. Especially useful in this context is the concept 
of weak differentiability [13]: the probability density p(X) is weakly differ- 
entiable if there exists a function s(-) G L 2 (p), such that for all / with 
/ \f(\)\ 2 p(\ + 6)d\ < oo the function g(9) = J f(X)p(X + 9)d\ has a deriva- 
tive g'(6) = J /(A)s(A + 6)p(\ + 9)d\. 

To get more insight into the structure of the optimal covariant observable 
for n > 1, let us consider it in the situation of the Example. Then 7i® n 
is isomorphic to L 2 (R n ), and A^ is just the operator of multiplication by 
A = Ai + . . . + A n . 

Theorem 2. The optimal covariant observable (16) is the spectral mea- 
sure of the selfadjoint operator 



Observable ^ (Q* <S> ■ ■ ■ <S> I + ■ ■ ■ + I <S> ■ ■ ■ <S> Q*) corresponds to a "semi- 
classical" method of estimation, when the optimal quantum estimates for 
each of n components in H® n are found and then used in a classical way to 
obtain the average over n observations. The theorem shows that projecting 
this average onto Ti.^ (and thus introducing entanglement) gives the optimal 
quantum estimate for n observations. 



g(n) = p wl (Q„ (g, ... (g, / + ... + / (g, ... (g, g,) pW. 



(18) 



n 
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Proof. As in the Example, one shows that (16) is orthogonal resolution 
of the identity. As follows from (17) the subspace = P^Ti. consists of 
the functions of the form 



<^(Ai, . . . , A r . 



c(A) 



V£W(A) 



^(Ai)...^(A„), 



(19) 



where j \c(X)\ 2 dX = \\(p\\ 2 . Consider the action on these functions of the 
self adjoint operator 



Q 



(n) 



xMi n \dx) 



I 



X- 



■ ixX \^ n \X)}(^ n \X')\e ixX ' dx 



||^)(A)||||^)(A0|| 2vr_ 
(Notice a parallel between this expression and the Pitman formula 

/ ep(x 1 -o)...p(x n -e)d9 



0* 



Jp(x 1 -9)...p(x n -9)d9 



(20) 



for the classical optimal covariant estimate of the shift parameter 9 in the 
family {p(x 1 — 9) . . .p(x n — 9)}), see e.g. [17]). By using the analog of (13), 
one shows that for ip e 



Qi n V(Ai,...,A n ) = i 



c'(A) 



V^W(A) 



^(Ai)...^(An), 



provided c(A) is absolutely continuous and c'(A) is square integrable. We 
have 

d_ 



i n ft 



hence 

gi n V(Ai,...,A n ) 



.y»n _d_ 

1 2^j=i ax. 



n 



c(Ai + . . . + A n ) 

v /p(n)( Al + ... + An ) 



^(Ai)...V(A„) 



-c(Ai + . . . + A n ) 



^(A!)...^(A n ) 
v/pW(A 1 + ... + A n ) 
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Taking into account that ip(\j) = ^/p(\j)e ia ^ Xj ^ and performing differentia- 
tion, we obtain after some transformations 

Ql"V(Ai,...,A n ) = V J -<p(Xi, ■ ■ ■ , Ki) 

n 



+-F(\ 1 ,...,X n )ip(X 1 , ...,X n ), 



where 



is a real function, or, briefly, 

f i 

Ql n V = - (Q* ® • • • ® / + • • • + / ® • • • ® <2*) y + T-fy, y? e 

n 2 
Taking inner product with ip, and noticing that both and 

f 

n 

are selfadjoint, we have (ip\Fip) — 0, ip e 7ii n ' ) , whence (18) follows. 

Let us also show directly that {(p\Fp) — 0, if </? is given by (19). We have 



(Q* <8> ... ® J + ... + J <g> ... <g> Q*) 



1 V W) (A)]' I^P'(^)\ |e(A)| 



<^> - y (^-^^i^.-.a)^,.^ (2 i) 

/ ,Jc(A)| 2 f f [p^(\)]' P (\ 1 )...p(\ n ) if 9 ^L-V 



But 



P (n) (A) = / p{\,)...p{\ n )d n - 1 a (22) 
7r(A) 

= p(A)* n = y'p(A 1 )...p(An-iMA-A 1 -...-A n _ 1 )dA 1 ...dA n _ 1 . 
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Differentiating the last equality with respect to A, we obtain 




[pW(A)]' = / p (Ai) . . .p(A n _i)p'(A - A! - ... - Xn^dX, . . . dA n _! 



Taking into account (22), (23) shows that the inner integral in (21) is equal 
to zero. □ 

4 The Limit Theorem 

In this Section we impose the regularity assumption: 

(A) ip e T>(A), or, equivalently, \\Aifj\\ 2 = J X 2 p(X)dX < oo, where p(X) = 
\\ip(X)\\ 2 is the probability density of observable A in the state S. Without loss 
of generality we assume E^(A) = (ip\Aip) = f Xp(X)dX = 0, then D S (A) = 
J X 2 p(X)dX is the variance of A. 

The uncertainty relation (5) together with (9) imply the lower bound for 
the variance of the optimal covariant estimate 



with equality attained if and only if p(A) is the Gaussian density The quan- 
tity D* (n) represents the accessible minimum of variances of arbitrary (in 
general, entangled) estimates of the shift parameter, based on n independent 
observations. 

By a variant of local central limit theorem [16], 



and similarly for all Xj. Hence 




(23) 




p n {X) := ^/nap(^fnaX) 



*n 



Po(A) : 



1 




V2n 
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in the sense of L 1 . Hence 



dX 



< J \ Po (X)-p n (X)\dX^O. 



Therefore the Fourier transform of y^p n (X) converges to that of >Jpo(X) in 
L 2 , and the probability density of the optimal covariant observable 

2 



Pg(n) (X) 



1 

2^ 



e- ixX y/p(\)* n d\ 



satisfies the local limit theorem 

r 2~ 



7T 



-2x z 



\fno 



P S ( n ) 



x 



y/na 



dx — > 0. 



Thus, under the assumption (A) the distribution of the optimal covariant 
observable in the state S$ is asymptotically normal with parameters 

(0, 4^2)- In this sense the bound of the uncertainty relation is asymptotically 
attainable. Moreover, by using the main result of [13] one has the asymptotic 
efficiency: 

Theorem 3. Ifp(X) is weakly differentiable and D* (n) is finite for some 
n, then 

nD *»> = iD^4)- 

5 Semiclassical Estimation 

Consider now the estimation strategy when the optimal covariant quantum 
estimate ilff ^ (dx) is found for every of the n components in the tensor prod- 
uct H® n , and then the classical estimation based on the obtained n outcomes 
is made. The probability density of observable ilff ^ (dx) in the state S d is 

1 2 

Pe(x) = —\ip(x-0)\ , 

where ip(x) = J e lXx \\ip(X)\\ dX. Under the assumption (A), ip(x) and hence 
pe(x) is differentiable, and for every unbiased estimate over n observations 
the classical Cramer-Rao inequality holds: 



D„ > 



n 



p e (x) 



dx 
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Notice that 



I %^x = 4 / l Re ^)y| 2 ^l = 4 / [mx)n * dx 
J Pe{x) J Wx)\ 2 2tt J Lin ;|J 2tt 



One the other hand, 



A 2 |h/>(A)|| 2 dA 



■0 (x) 



2 dx 
2tt' 



Comparing this with the asymptotically attainable quantum bound of the 
uncertainty relation, we have 



/ 



■0 (x) 



2 <ix 
2^ 



2 dx 
27' 



if only f3(x) := arg-0(x) ^ const. Thus, under this condition, for n large 
enough 

nmin D n > nD* 



where the minimum is over all unbiased classical estimates using unentan- 
gled quantum observables, which demonstrates superiority of the entangled 
quantum estimation. 



6 Irregular Case 



We now consider an instance of the Example where the regularity assumption 
(A) does not hold. Let S = with 

1 sin aA 

^na A 

so that Dg(A) = oo. In the coordinate representation this corresponds to the 
rectangular function 

fa) = [ e-^{\)d\ = i V^' .'f X f. [_a ' a]; 
J I 0, ifx ^ [-a, a], 
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that is, to the particle position uniformly distributed in [—a, a], with the 
probability density p(x) = ^\ip(x) | 2 in R. Unfortunately the probability 
density of the optimal covariant observable Q* which is 



2 

1 J e- iXx \i){\)\d\ 



cannot be found in explicit form, therefore we shall consider semiclassical 
estimation based on unmodified position observable Q having the probability 
density p(x — 9). 

Turning to the case of n observations, we denote 

Qi = Q <8> ... ® Q n = I <E> ■■■ ® Q- 

As is well known (see e. g. [2], n. 28.6), the Pitman estimate (20) in the 
case of the rectangular probability density p(x) has the form 

9* = ^ [min (Q 1 , Q n ) + max (Q u . . . , Q n )\ . 
Its variance is equal to 

2a 2 2a 2 



Ds (0* 



(n + l)(n + 2) n 



2 



which shows faster decay than 1/n characteristic to the regular case. The 
rectangular density is nondifferentiable violating the regularity assumption 
and making possible more efficient estimation than one which would follow 
from the Cramer- Rao bound (inapplicable in this case). Moreover, denoting 
by p n (x) the probability density of 9*, one has the limit law [2] 

j \x | 

lim n~ l p n (x/n) = —e~~. (24) 

rwoo 2d 

Now we shall find the asymptotics of the optimal quantum covariant 
estimate (18). Denoting by 

/ . , ,/>\,2 1 /sinaA x z 
p(A) = |^(A)| 2 = — ( — — 
ixa \ A 
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the probability density of A, and by f(x) = f e tXx p(\)d\, we have 



/(*) = 



fl-g, ifxe[-2a,2a]; 



Thus 

whence 

lim np n (n\y 



0, if x ^ [—a, a]. 
1 



p( A )« = - / e»*f(xrdx, 



— lim / e a */(:r/n) n d:r 

27T n->oo J 

— /exp | iAx — | dx = - — — 7— — r — — . 

2vr y F V 2a / 27ra A 2 + (2a)- 2 



Taking into account that 

~ iXx - 1 =d\ = 2K f M . . 

where i\~o is the Macdonald (modified Bessel) function, and by using the limit 
law, we can show (see Appendix) that the renormalized probability density 
of the optimal observable 



Plin) (X) - — 

obeys the limit law 



e- iXx y/p(\)* n d\ 



2 

na 



^ ( f- 

2a 



dx = 0. 



(25) 



With some more effort we can show (see Appendix) that its variance satisfies 



lim n 2 D* s(n) = lim / ( — \J np{n\)* n J d\ 

n^oo n^oo \d\ J 



(26) 



d 



2ira J K \d\ U A 2 + (2a)" 



dA 



2 ' 
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which is four times less than for the semiclassical estimate 6*. Moreover, 
K (x) ~ \fj^ e ~ x f° r large positive x, whence 

2 _ n 
\x\ 

which shows that the tails of the asymptotic distribution of the optimal 
estimate has somewhat faster decay than (24). 
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Proof of (26). Denoting f (x) = e~W 2a , p (X) = 1 



2na X 2 + (2a)-' 2 ' 

Mn{x) = fo(x) - f(x/n) n ; Ap n (\) = p (X) - np(n\y 



we have 



Ap, 



e iXx Af n (x)dx. 



Now we observe that 



lim / \Af n (xY k) \ \x\ l dx = 0; k,l = 0,1,... 



From (27) it follows in particular that 



max|Ap n (A)| < — 

A Z7T 



| Af n (x)\ dx — > 0, asn — > oo. 



(27) 



(28) 



Let us also estimate the tails of Ap n (X). Taking into account that A/ n (0) 
(A/ n (0))' = 0, and making twice integration by parts in (27), we obtain 



/ 


A/„(x)" 







(29) 



where lim^oo e n = 0. In the same way we obtain 



max 

A 



< J \xAf n (x)\ dx — > 0, asn^oo, (30) 
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and 



^A P „(A) 



< 7VT3 , 
|A| A^o 



|(xA/ B (x)r|da; 



|A 



3 ' 



(31) 



From (28), (29) it follows that J |Ap n (A)| d\ — > as n — > oo, hence, arguing 
as before theorem 3, we have ^/np(n\)* n — > a/po(A) in L 2 , implying (25 ). 
The estimates (28), (29), (30), (31) together with 



Po(A) 



1 + A 2 ' 



dA 



Po(A) 



1 + IAI 



imply 



d \ ^ ( d ^ ^ 



K(A) 2 K(A) 2 



Po(A) p n (A) 



e" 



1 + A 2 ' 



whence (26) follows. 
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